Controllable Coupling between Flux Qubit 
and Nanomechanical Resonator by Magnetic Field 
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We propose an active mechanism for coupling the quantized mode of a nanomechanical resonator 
to the persistent current in the loop of superconducting Josephson junction (or phase slip) flux qubit. 
This coupling is independently controlled by an external coupling magnetic field. The whole system 
forms a novel solid-state cavity QED architecture in strong coupling limit. This architecture can be 
used to demonstrate quantum optics phenomena and coherently manipulate the qubit for quantum 
information processing. The coupling mechanism is applicable for more generalized situations where 
the superconducting Josephson junction system is a multi-level system. We also address the practical 
issues concerning experimental realization. 

PACS numbers: 74.78.Fk, 85.85.+J, 03.65.Lx 



I. INTRODUCTION 

In recent years, great advances in improving the co- 
herence of superconducting qubit have made it a promis- 
ing candidate for the physical realization of quantum in- 
formation processing. Single qubit Rabi oscillation and 
Ramsey fringe have been observed and two qubits en- 
tanglement is also achieved. Meanwhile, as an artificial 
two-level atom, superconducting qubit is adjustable (e.g. 
by flux, bias voltage, etc.) and scalable. These features 
are favorable for quantum state engineering. A number 
of protocols are proposed to engineer the superconduct- 
ing qubit to form a quantum network. Among them, a 
very intriguing and successful example is the circuit QED 
architecture [l| . By coupling the Cooper pair box (charge 
qubit) to the quantized field of a coplanar superconduct- 
ing transmission line, a macroscopic solid-state analog of 
cavity QED is realized on chip. Most recently, vacuum 
Rabi oscillations are observed in a coupling system of 
3- J J flux qubit and LC circuit @ . Quantum optical phe- 
nomena in traditional cavity QED can be demonstrated 
in this solid state composite system. Furthermore, due 
to its special structure, it offers a number of advantages, 
such as strong coupling and easy controllability. Thus, 
some protocols that cannot be realized previously in the 
optical cavity QED now becomes possible [1, 0| . 

The circuit QED experiments motivate us to investi- 
gate the possibility of substituting other quantum solid 
state devices for the transmission line. It is much desir- 
able to couple Josephson junction qubit to a device with 
low energy consuming and small size. If the strong cou- 
pling and easy controllability can also be achieved, we get 
another favorable cavity QED structure. A possible can- 
didate of this solid state device is the nanomechanical res- 
onator (NAMR). Nanomechanical resonators of GHz os- 
cillation have already been observed. It is supposed that 
the nanomechanical resonator enter the quantum regime 
at the attainable temperature of the dilution refrigerator. 



The schemes of coupling Josephson charge qubit or phase 
qubit to NAMR have already been proposed. Based on 
these coupling mechanisms, several quantum state engi- 
neering protocols are put forward [ESjHJs, 9]. However, 
due to the difficulty to reach quantum regime of NAMR, 
those protocols have not been implemented experimen- 
tally yet. On the other hand, the coupling mechanism of 
NAMR and flux qubit is also an attractive problem since 
the flux qubit is supposed to have longer coherence time 
as it is less affected by the charge fluctuation in the struc- 
ture. To our best knowledge, this has not been studied in 
detail previously. Here, we present a novel mechanism of 
coupling NAMR to flux qubit. As we present below, the 
coupling strength between the NAMR and the flux qubit 
can be adjusted conveniently by a coupling magnetic field 
and turned on and off within the coherence time of flux 
qubit. Since the coupling magnetic field is independent of 
the single qubit operation, it is possible to make it strong 
enough even for GHz oscillation. Therefore our proposal 
is in principle possible to approach the "strong coupling 
regime" of cavity QED at attainable temperature of dilu- 
tion refrigerator. This coupling system acts as an analog 
of cavity QED system with more flexibility. We expect 
that it enables various applications to quantum informa- 
tion processing and quantum state engineering. 

This paper is organized as follows: in SecHIl we briefly 
review the setup of 3 Josephson-junction flux qubit and 
NAMR as well as their experimental progresses. Then 
we get into the coupling mechanism for flux qubit and 
NAMR. This coupling mechanism can be equally ap- 
plied to rfSQUID flux qubit and phase slip flux qubit. In 
Sec JIIH the spectrum of the coupling system is presented 
in the "weak coupling" and the "strong coupling" limits 
respectively. The readout and quantum nondemolition 
measurement for the flux qubit is studied in Sec lIVl Wc 
also consider the application of this coupling mechanism 
in quantum computation in Sec|V]and generalize the cou- 
pling system to beyond spin-boson model in Sec lVIl The 
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possible problems on experimental realization and their 
solutions are given in Sec. IVIII In Sec lVIIIl some discus- 
sions and remarks are included. 



II. COUPLING NANOMECHANICAL 
RESONATOR AND FLUX QUBIT BY A 
TRANSVERSE MAGNETIC FIELD 

A. The 3-junction flux qubit 

Josephson charge qubit system has been used to cou- 
ple with NAMR. Here, we study another superconduct- 
ing qubit system - 3 Josephson junction (3-JJ) flux qubit 
[ToL ITD . EU EH- I n contrast with charge qubits, the 
flux qubit is far less sensitive to charge fluctuations. Es- 
timations show that the flux qubit have a relatively high 
quality factor [ll|, [T^]. The configuration of flux qubit 
consists of a superconducting loop with three Josephson 
junctions, and the Josephson coupling energy is much 
larger than the charging energy for each junction. The 
quantum state of this system is mainly determined by 
the phase degree of freedom. The Josephson energy of 
the three Josephson junction loop reads 

U(ipi,ip 2 )= - Ej cos <px - Ej cos ip 2 

- aEj cos (2nf - y x - Vi) , (1) 

where the constraint of fluxoid quantization has already 
been taken into account. Here, Ej is the Josephson 
coupling energy of two identical junction and tpi, ip 2 
are phase differences across the two junctions respec- 
tively. The Josephson energy of the third junction is 
aEj , f — /$o with $ f the external flux applied in the 
loop and <£>o = h/2e the flux quantum. In the vicinity 
of / = 0.5, if a > 0.5, a double-well potential is formed 
within each 2ir x 2ir cell in the phase plane and the two 
lowest stable classical states have persistent circulating 



currents I p = 2eEjy 1 — (l/2a) /h with opposite direc- 
tions. Therefore, the flux qubit is also called persistent 
current qubit. Within the qubit subspace spanned by 
{|0) , |1)} (|0) and |1) denote clockwise and counterclock- 
wise circulating states respectively), the Hamiltonian of 
the qubit system reads as 



Hf = ujf<7 z + Aa x = Q,5 Z , 



(2) 



with u>f = Ip&o (/ — 0.5) is the energy spacing of the 
two classical stable states and A the tunneling split- 
ting between the two states, fl — \Ju"j + A 2 and a z = 

cos 8a z + sin 9a x , tan9 = A/u>f. The offset of / from 
0.5 determines the level splitting of the two states and 
the barrier for quantum tunneling between the states 
strongly depends on the value of a. If the third junction 
is replaced by a dc SQUID, both / and a are tunable in 
experiments by the applied flux or the microwave current 

Halm. 



B. The nanomechanical resonator 

The flexural modes of thin beams can be described by 
the so-called Euler-Bernoulli equations [TBI ]. In our pro- 
posal only the fundamental flexural mode of the NAMR 
is taken into account. All the other modes have a 
much smaller coupling to flux qubit and can be ne- 
glected [l(| E3] ■ In this case, the NAMR is modeled as a 
harmonic oscillator with a high-Q mode of frequency ojb- 
The Hamiltonian without dissipation reads @, 0| 
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H = £*. + -mu> 
2m 2 



(3) 



In pursuing the quantum behavior of macro scale object 
the nano scale mechanical resonator plays an important 
role. At sufficient low temperature the zero-point fluc- 
tuation of nano mechanical resonator will be comparable 
to its thermal Brownian motion. The detection of zero- 
point fluctuation of the nano mechanical resonator can 
give a direct test of the Heisenberg's uncertainty prin- 
ciple. With a sensitivity up to 10 times the amplitude 
of the zero-point fluctuation, LaHaye et al. have experi- 
mentally detected the vibrations of a 20-MHz mechanical 
beam of tens micrometers size [l8[. For a 20-MHz me- 
chanical resonator its temperature must be cooled below 
1 mK to suppress the thermal fluctuation. For a GHz 
mechanical resonator a temperature of 50 mK is suffi- 
cient to effectively freeze out its thermal fluctuation and 
let it enter quantum regime. This temperature is already 
attainable in the dilution refrigerator. 

The lithographic technology for NAMR is rather ma- 
ture. The important advantages of NAMR are the poten- 
tially higher quality factor and frequency comparable to 
superconducting qubit. Ever since the early demonstra- 
tion of a radio frequency mechanical resonator at Caltech 
[19] . great advances have been made. The attainable fre- 
quencies for the fundamental flexural modes can reach 
590 MHz for the doubly-clamped SiC mechanical res- 
onator of the size 1 x 0.05 x 0.05 /im [2(| and 1 GHz 
oscillation frequency has also been measured [2lj . It is 
argued that quantized displacements of the mechanical 
resonator were observed despite of some opposite opin- 
ions |22l |. For a 1 /im beam a quality factor Q of 1700 
has been observed at a frequency of 110 MHz [23j |. In a 
carefully designed antenna shape, Gaidarzhy et al. have 
achieved Q — 11000 for 21 MHz oscillation at the temper- 
ature of 60 mK and Q = 150 for 1.49 GHz oscillation at 
the temperature of IK with a comparatively large double 
clamped beam (24| . The significantly small size (~ /im) 
of the NAMR is also favorable for incorporating it in the 
superconducting qubit circuit. 



C. The composite system with tunable coupling 

To achieve a "strong" interaction, the coupling dynam- 
ical variable usually should be the dominant one in the 
dynamics of the composite system. For the Josephson 
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FIG. 1: (Color on line) A 3-Josephson-junction flux qubit loop 
is located in the x-y plane and a NAMR is integrated in the 
loop(indicated by a green line). The z-direction oscillation 
of the NAMR couples to the current in the flux qubit loop 
by a transverse magnetic field Bo in y-direction. Another 
tunable magnetic flux penetrates this loop tunes the free 
Hamiltonian of the 3-JJ system. 



phase qubit [25|, |26| , the phase degree of freedom domi- 
nates the dynamics and the bias current coupled with the 
phase is modified by the dilatational motion of the piezo- 
electric dilatational resonator [27[ . While for the Joseph- 
son charge qubit, the coupling mechanism is that the 
resonator displacement modifies the effective bias charge 
of a Cooper-pair box [3, H, [13, HH . These previous inves- 
tigations enlighten us to consider the coupling between 
persistent current in superconducting flux qubit loop and 
the motion of nano mechanical resonator. 

Since the Josephson coupling energy of each junction 
in the flux qubit is much larger than that in the charge 
qubit, the persistent current in the loop could be about 
hundreds of nano ampere fl2| in contrast with the critical 
current of the charge qubit (usually about 20 ~ 50 nA). 
The magnitude of this persistent current naturally leads 
us to consider the magnetomotive displacement actua- 
tion and sensing technique pHHH. It is well known when 
a current passes through a beam with conducting mate- 
rial, the perpendicular arrangeed of an external magnetic 
field and the direction of the current generates a Lorentz 
force in the plane of the beam. This is just the actuation 
part of the magnetomotive technique. Meanwhile, the re- 
sulted displacement of the beam under the Lorentz force 
generates an electromotive force, or voltage, which serves 
as measurement. Thus, if the doubly-clamped nano beam 
coated with superconducting material is incorporated in 
the superconducting qubit loop, the persistent current in- 
duces a Lorentz force with opposite directions for clock- 
wise and counterclockwise current. The oscillation of the 
NAMR is modulated by these Lorentz force. In this way, 
the quantized harmonic oscillation mode of the beam is 
coupled to the quantum state of the flux qubit system. 
And this is just the coupling mechanism considered in 
our paper. 

Our proposal is illustrated in Fig[T] A 3-JJ system is 
fabricated on the x-y plane. The external applied mag- 
netic flux is enclosed in the loop modulated by the 



control lines (the lines are not plotted). The 4- J J ver- 
sion of flux qubit system can also be used here to allow 
the modulation of the effective Josephson energy of the 
third junction and hence the tunneling amplitude of the 
two current states. One side of the loop (indicated by 
thick (green) rod) is suspended from the substrate and 
clamped at both ends. This can be fabricated with a 
doubly-clamped nanomechanical beam coated with su- 
perconductor or with the superconductor itself as the me- 
chanical resonator. A magnetic field Bo is applied in the 
y-direction. As we discussed above, the circulating sup- 
percurrent under the magnetic field generates a Lorentz 
force in the z-direction. The magnitude of the force is 
BqI p L, with L the effective length of the resonator along 
x-direction (L = £Lq and Lq is the actual length of the 
resonator, £ a factor depending on the oscillation mode 
[29l | . for the fundamental oscillation mode of a doubly 
clamped beam £ « 0.8). This force results a forced term 
in the Hamiltonian, which reads Hfb = Fz — BqI p Lz. 
With the two-level approximation of 3-JJ loop and the 
singlemode boson approximation, the coupling is written 
as 



H fb = g(a + a))cr z , 
for z ~ a + . Here, 

g = B {t)I p LS z 



(4) 



(5) 



and S z = y/h/2mu>b is the amplitude of zero point motion 
in z-direction of the NAMR, with m the effective mass 
of the resonator, Ub the frequency of the fundamental 
nexural mode; a (af) is the creation (annihilation) op- 
erator of the mode of the flexural motion in z-direction. 
0z = |0) (0| — |1) (1| is the Pauli matrix defined in the 
basis of |J0) , |1)}. We see that this interaction Hfb ac- 
tually couples the two systems. Together with the free 
Hamiltonian of flux qubit and NAMR, the Hamiltonian 
of whole system reads 



H = u>bd ;'<z + ujfcr z + Aa x + g [a + a') a 



(6) 



An important advantage of this coupling mechanism 
is the convenient controllability. As seen from Eq.©, 
the coupling constant is directly dependent on the ap- 
plied coupling magnetic field B . Thus, both the magni- 
tude and sign of the coupling constant can be modified. 
What's more important is that the control parameter Bq 
in coupling coefficient Eq. ^ is independent from the pa- 
rameters of free Hamiltonian, such as bias voltage and 
external magnetic flux $ f . This means the free Hamilto- 
nian and the interaction Hamiltonian can be manipulated 
independently. This full controllability is a rather favor- 
able feature for quantum state engineering and quantum 
information processing protocols. It is in contrast with 
the coupling of charge qubit and NAMR, where the cou- 
pling strength is controlled by the bias voltage which is 
also the crucial parameter to determine the energy spac- 
ing of the charge qubit. For example, for bang-bang 
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cooling of NAMR by charge qubit |8|, the bias voltage 
should be set to certain value to induce desirable damp- 
ing. Therefore the on-and-off of the interaction between 
the qubit and NAMR can only approximately controlled 
by detuning and this can result in harmful reheating of 
the NAMR. But in our present coupling mechanism, both 
the coupling coefficient and the energy spacing are inde- 
pendent. Thus, this "bang-bang" cooling protocol should 
be implemented more reliable by the flux qubit and the 
NAMR with above coupling mechanism. 

To estimate the coupling strength, we use the follow- 
ing parameters in refs. @, Ho, [3lj] : Ip = 660 nA, L = 3.9 
fim, uj b = 100 MHz, 6 Z = 2.6 x 10" 13 m, Q = 2 x 10 4 
and assume the applied magnetic field to be Bq = 5 mT. 
Then we have g sw 4.01 MHz. Hence we see that the 
"strong coupling" regime for cavity QED is potentially 
realizable in our scheme. This regime requires the period 
of the Rabi oscillation l/g, is much shorter than both 
the decoherence time 1/7 of the two- level system and 
the average lifetime 1/k — Q/u>b of the "photon" in the 
"cavity" 32j . For this composite system, the decoherence 
time for flux qubit is 1 — 10 /is and the cavity lifetime is 
about 200 /its, while the Rabi oscillation time 0.016 fis is 
much shorter than the two lifetime scale. For GHz oscil- 
lation, the quality factor is rather low[24| (1.49 GHz with 
Q = 150). This corresponding to a much shorter cavity 
lifetime (about 0.1 /is. However, the coupling strength 
can be increased by larger coupling magnetic field. For 
example, if we take Bq = 50 mT, the Rabi oscillation pe- 
riod l/g w 0.016 /Lis which is still short enough to reach 
"strong coupling regime" . Therefore, this protocol might 
be promising in dilution refrigerator (several tens of mil- 
likelvin). 



D. Phase slip flux qubit and NAMR 

Most recently, a new type of flux qubit - phase slip 
flux qubit is proposed based on coherent quantum phase 
slip [33|, Hi]]. Phase slip flux qubit is formed by a 
high-resistance superconducting thin wire instead of the 
Josephson junctions. The computational basis are also 
the two opposed persistent current states. Our coupling 
scheme can be equally applicable to this type of flux 
qubit. There are two important advantages to consider 




FIG. 2: (Color on line) A superconducting thin wire in the 
loop acts as the center of phase slip flux qubit and a NAMR. 



phase slip flux qubit. First, since the superconducting 
thin wire can be fabricated by a suspended carbon nan- 
otube, it acts as the phase slip center and the NAMR si- 
multaneously. The circuit configuration is simplified (see 
Fig. [2]). Secondly, if one use this qubit, one can free from 
any fluctuator due to imperfection or two-level systems 
hidden in the dielectric layer of Josephson junction. 



III. ENERGY SPECTRUM OF COUPLING 
SYSTEM 



The Larmor frequency of superconducting qubit is 
about the order of 10 GHz, while the frequency of NAMR 
only reaches several hundred MHz with quality factor 
10 4 at present stage. Thus, the composite system of flux 
qubit and NAMR is in large detuning regime of cavity 
QED, i.e., the following condition is satisfied 



9 



\n - 



< 1. 



(J) 



However, the superconducting flux qubit and the NAMR 
is non-resonant, i.e. Q LOb- This is in contrast with 
Yale's circuit QED experiment where the C oop er pair 
box is resonant with the ID transmission line [351 ]. In the 
following, we discuss the energy spectrum of our model 
in two different regimes: g -C u>b (denoted as "weak cou- 
pling") and g w LUb (denoted as "strong coupling"). In 
our proposal, the two regimes can be reached by varying 
the applied coupling magnetic field Bq. And the energy 
spectrum are qualitatively different from each other. It 
is notable that Ref.[36[ has proved that the dispersive 
measurement back action can be enhanced or reduced by 
cavity damping respectively in the two regimes. 



A. "Weak coupling" and sideband spectrum 

For the parameters B = 5 mT, g = 4.01 MHz, both 
g/n and g/cub are much smaller than 1. In this case, the 
energy spectrum can be calculated by Floquet approach 
or by Frolich transformation. After performing a unitary 
transformation on the original Hamiltonian ([5]), we get 
the effective Hamiltonian 



H e ffi ~ ujba) a + Qa z + i 



.g 2 sm26 



2 _ „t2\ X 



g sin ( 
0~~ 



(a + a^) 2 (cos9a x + sm9a z ) . (8) 



The spectrum are nu>b + mfl plus some small off-diagonal 
transition terms that are of order 0(g/u>) or 0(g/il). 
As shown in FigfSl the energy levels of the two subsys- 
tems are weakly perturbed by the coupling due to the 
large detuning. By applying microwave pulse to induce 
the transition between those levels, the blue sideband 
(| 00) -> |11)) and red sideband transition (|01) -> 1 10>) 
can be observed in addition to the main zero-photon tran- 
sition 1 00) — > 1 10) [13|- This atomic physics phenomenon 
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FIG. 3: (Color on line) The dressed energy level and transition 
diagram for the weakly coupled system. 

has already been observed in solid quan tum system of 
flux qubit and dc SQUID oscillator [33]. For our pro- 
posal, similar spectra are expected. 



B. "Strong coupling" and dispersive shift 

With large magnetic field, for example, Bq = 100 mT, 
then g = 80.2 MHz. The magnitude of coupling is com- 
parable to the characteristic energy scale of the NAMR. 
In this case, the coupling term is not a perturbation with 
respect to the free Hamiltonian of NAMR. Therefore, we 
can not use perturbation theory. However, since the large 
detuning condition is still hold, we resort to adiabatic 
elimination (or coarse-graining technique) to deal with 
this problem. Since £1 3> the energy spectrum of the 
whole system is roughly energy band structure. Then the 
spacing of bands is determined primarily by and the 
energy spacing within each band is approximately u>b- 
By adiabatic eliminating the transition between differ- 
ent bands, we can obtain the effective Hamiltonian from 
eq.© 



H = H e \e) (e\+H g \g) (g\, 

where the two component Hamiltonians are 

2 a <fb_ 

°lg 



H e .. 



jAl A^g ±n-g 2 cos 2 9 



with the frequencies 



2 ig 2 uJb sin 2 



J ll — UJb 

and the new bosonic operators A e g are defined by 
In the above equation 



/J>e,g 
Ve,g 
Ve 



/uJ e,g 

4w fc 



OJb 



e,9 



4uj, 



e,9 



g cos 6y/ujbUJ e 

9 ' 

g cos 9 y/u>bUJg 



(9) 
(10) 

(11) 
(12) 

(13a) 
(13b) 
(13c) 
(13d) 



i 
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FIG. 4: Dispersive "pull" of the frequency of NAMR with 
asymmetric shifts from its original frequency. 



Thus, for different state of the flux qubit, the energy 
spectrum of the NAMR is shifted by different value, or 
in other word, the flux qubit pulls the cavity frequency 
by Si and —82 (see FigQ| with 



Si 



4g 2 ujb sin 2 9 



Q — ujb 



■ u>b, 



j 2 u>b sin 2 1 

— UJb 



(14a) 



(14b) 



In contrast with the dispersive limit discussed in circuit 
QED [l[, the dispersive shift here is asymmetric. 



IV. QND MEASUREMENT FOR FLUX QUBIT 

Usually, the flux qubit is measured through a dc 
SQUID in the underdamped regime that is inductively 
coupled [U,[39| or directly coupled 0, lU, E2 . The ex- 
ternal applied flux plus the flux induced by the persistent 
current decide the switching current of the dc SQUID. By 
ramping a bias current to the dc SQUID, the switching 
current can be recorded. When a continuous microwave 
is resonant with the level spacing of the two eigenstates 
of the flux qubit, the qubit is flipped and the switch- 
ing current is changed. This results in peak and dip in 
the switching current level versus the applied flux. With 
this method, both the energy spectrum and the dynamic 
evolution have been observed. During this measurement 
process, the quantum information encoded in the qubit is 
destroyed (30|. It would be favorable to design a nonde- 
structive and quantum nondemolition (QND) [431 ] mea- 
surement protocol. A novel phase-sensitive microwave 
reflection approach are no w appli ed for the readout of su- 
perconducting qubits [13, 144L l45l . H(| . The advantage of 
this method is that it directly probes the dynamics of the 
Josephson plasma resonance in both the linear and non- 
linear regime without switching the detector Josephson 
circuit to dissipative state. It succeeded to provide very 
fast and far less destructive measurement of the qubit. 
However, the QND readout for the flux qubit [47[ only 
works far away from the optimal point, where the qubit 
coherence is destroyed very quickly. 
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Here, we indicate that our composite system can be 
used to perform QND measurement on flux qubit. As we 
discussed in the previous section, in both cases, the inter- 
action between the NAMR and the flux qubit results in 
the mixed energy spectrum for them. Therefore, through 
the spectroscopy measurement, the quantum state of one 
system can be detected via the spectrum of the other one. 
Especially, if the interaction Hamiltonian commutes with 
the free Hamiltonian of the measured object and does not 
commute with that of the measuring device, a QND mea- 
surement protocol can be implemented. This is just the 
case of "strong coupling" limit in our proposal. In this 
case, by the spectroscopy of the NAMR, the flux qubit 
state can be read out without perturbation. This can 
be predicted from Eq.([9HT2"]): the frequency of NAMR 
is u> e when the flux qubit is in the excited state, while 
its frequency is u g when the flux qubit is in the ground 
state. The interaction between the NAMR and the flux 
qubit commutes with the free Hamiltonian of the flux 
qubit. Hence the measured probabilities of eigenstates 
are not perturbed by this readout. Therefore, the spec- 
troscopy measurement of NAMR provides a high resolu- 
tion QND measurement on the qubit state. It should be 
noticed that this scheme does not work at the exact op- 
timal point (sin 9 = at the optimal point). However, if 
only the operation point is a little bit shifted away from 
the optimal point, e.g. u>f rs A, the resulting frequency 
shift is observable (suppose Bq — 100 mT, u>b = 100 MHz 
and Q = 2 x 10 4 , then Sip ^> k) in the spectrum. 

This frequency shift can be measured through the fre- 
quency measurement of the NAMR. In pri ncip le, this can 
be done with magnetomotive technique [2J, [48|. Dur- 
ing the measurement process, a perpendicular magnetic 
field and oscillation current is applied on the NAMR. 
Then, the NAMR behaves like a frequency dependent re- 
sistance. The largest effective resistance is obtained when 
the NAMR is resonance with the oscillation current. 
Thus the frequency of NAMR is inferred by the resonance 
peak of its voltage between its two end when we vary the 
frequency of the oscillation current in it. Another possi- 
ble way to the frequency measurement of NAMR is to use 
a single-electron transistors (SET) [H US [HI . The SET 
does not require extra magnetic field which might induce 
unwanted perturbation to the superconducting loop. It 
is supposed to have very high sensitivity and is expected 
to reach the limit by the uncertainty principle. For this 
method, the mechanical motion of the NAMR couples to 
the SET through a lead on the NAMR that is close to 
the island of SET. The motion of the NAMR modulates 
the coupling capacitance between the lead and the SET. 
When there is a bias voltage on the lead of NAMR the po- 
tential of the island near the NAMR is modulated. The 
frequency of the mechanical motion is detected through 
the conductance of the SET. 

On the other hand, the measurement on high- 
frequency mechanical oscillator has some practical diffi- 
culties as additional strong coupled transducer is required 
to convert its dynamics to electronic signals. There is an 
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FIG. 5: (Color on line) The NAMR couples to (a) two or (b) 
multiple flux qubits. The green narrow box represents the 
NAMR and the cross stands for the Josephson junction. The 
coupling is controlled by a magnetic field perpendicular to the 
NAMR in the coplanar of the NAMR and flux qubit. 



equally intriguing problem: to detect the property of the 
NAMR via the measurement of superconducting qubit 
since there have been some good measurement protocols 
of the latter. For the coupling mechanism presented in 
this paper, the effective Josephson energy is modified by 
the displacement of NAMR in a similar way as the charge 
energy of SET being modified by the displacement of 
NAMR. Thus, the flux qubit might be able to act as a 
transducer to detect the state of NAMR. Most recently, 
a QND measurement for NAMR via rf SQUID has been 
considered based on a configuration similar to ours [52| . 



APPLICATIONS IN QUANTUM 
COMPUTATION 



One of the possible applications of our proposal in 
quantum computation is to couple two or more flux 
qubits together and to realize two qubit logic gate. As 
shown in FigJS] (a) , the NAMR serves as a quantum data 
bus and the two identical qubit loops are connected to 
it. The total Hamiltonian reads 



H = uj b a^a + V" {ujj<7 zi + Aa xi 



a' er. 



(15) 



where a Z i and Pauli operators for the i th flux 

qubit. The coupling coefficient g can be modulated by 
the magnetic field Bq. If we fix Bq, the parameters for 
manipulating single qubit (i.e., ujf and A) operations are 
still tunable by adjusting / and a. This offers a universal 
architecture to realize coherent two qubit quantum logic 
gate [Hj]. Considering the inductance of the loop L p , 
we note that there is also direct magnetic coupling in- 
duced by the sharing edge L p (Ii +I 2 ) 2 /2 ~ a z ia z2 , here 
Ji, Ii are the currents in the 1st and 2nd qubit loop re- 
spectively. The order of magnitude of this zz coupling is 
about 10 MHz ~ 100 MHz. At the degeneracy point, this 
always-on coupling is commute with the whole Hamilto- 
nian, therefore it is not very hard to deal with. Away 
from the degeneracy point, the direct magnetic coupling 



7 



may have positive or negative effect with respect to con- 
crete proposals. 

More qubits can also be connected in the same way as 
shown in Figf5] (b) . By the tunable energy spacing of the 
flux qubit, we can selectively couple two qubits through 
the NAMR. With this configuration, we can realize the 
logic gate of two arbitrary flux qubits by adiabatically 
elimination o r by d y namically cancelation of the NAMR 
cavity mode [Mill. 



VI. BEYOND SPIN-BOSON MODEL 

In the above discussion, we have described our cou- 
pling mechanism in spin-boson regime where the 3-JJ 
superconducting loop is treated as a quasi two-level sys- 
tem, i.e., a qubit. However, when the magnetic flux is 
tuned away from $0/2, the lowest two energy levels can- 
not be isolated from other energy levels (see the energy 
spectrum, for example, in ref.[55||. Taking more energy 
levels into consideration is advantageous to investigate 
many intriguing phenomena that are traditionally stud- 
ied in atomic physics and quantum optics. For example, 
with the lowest three energy levels, stimulated Raman 
adiabatic passage (STIRAP) can be studied and some 
interesting behaviors, such as electromagnetic induced 
transparency (EIT) and dark states, can be exhibited 
in the coupling system. What's more, as the symmetry 
and the selection rule of the 3-level superconducting loop 
are different from that of the 3-level natural atom, some 
novel features can be demonstrated, for example, A-type 
atom [56| and persistent single photon generation [57j . 

Our coupling protocol can also be generalized to the 
model of multi-level atom in cavity. In this case, the 
dynamics of the 3-JJ loop is not confined to the two- 
level subspace. With newly-defined variables ip p — 
{ipx + Lpi) /2, ip m — (ipi — Lp 2 ) /2 and their conjugate mo- 
mentum P p , P m , the free Hamiltonian of the 3-JJ system 
has a similar form to that of a particle in a 2D periodical 
potential [Tlj : 

H f = i + ^+2Ml-co S ^co^ m ) 

+aEj (1- cos (2nf + 2<p m )), (16) 

where M p = 2Cj ($ /2tt) 2 and M m = M p (1 + 2a), and 
Cj is the capacitance of the first and second junctions. 
Then the interaction Hamiltonian is still Hfi, = BqILz 
where I is the current flows through the NAMR [57| 

1 = 6 S 3 ( 2 cos sm - sin (27r/ + 2ip m )) . (17) 
raw 

with 1/Cs = J2i l/Cj- The current I induce the transi- 
tion between different eigenstates of the free Hamiltonian 
of 3-JJ system and numerical calculation can predict the 
transition amplitudes. For example, if we only consider 
the lowest three energy levels |0), |1) and |2), the above 



0.4 




0.0-1 1 , 1 , , , 1 , 1 , 1 , ,— 

0.47 0.48 0.49 0.50 0.51 0.52 0.53 



FIG. 6: (Color on line) The transition amplitudes between 
the lowest three energy levels |0), |1) and |2) vary with the 
offset /. 



Hamiltonian can be written in a three dimensional sub- 
space: 

H = WbdW^Oj \i) (i| + A^(fiij K) 01 + h.c) (a + a r ) 

i i=jLj 

(18) 

where = 0,1,2) is the i-th eigen level of the 3-JJ 
system and f2; is the corresponding eigenenergy (usually 
we take f2 = 0). The coupling coefficient A = B (t)L6 z 
and flij = (i\I\j). As shown in FigJBl the transition 
amplitudes f^oi, ^02 and Q12 between the lowest three 
energy levels depend on /. That is to say, we can con- 
trol the transition between different energy levels of the 
uncoupled system. Apparently, this feature can be used 
in STIRAP technology and some other pro pos als based 
on three- level atom with the quantized field [58j . For ex- 
ample, when the superconducting loop is biased a little 
bit away from the the optimal point, a A type three- level 
atom is formed approximately. Following the same strat- 
egy of ref. [5^| , replacing the lowest three energy levels of 
quantronium by the eigen-levels considered above, Fock- 
states of the NAMR can also be generated. 



VII. EXPERIMENTAL CONSIDERATIONS 

In the above discussion, we have assumed that the 
NAMR oscillates only in the z-direction. However, to 
manipulate the 3-JJ flux qubit, a bias magnetic flux 
is applied through the loop in the z-direction. This bias 
magnetic field also induces a Lorentz force on the NAMR 
in y-direction. Therefore there exists an always-on inter- 
action between the flux qubit and the y-direction motion 
of the NAMR H'^ b = g' (a v + aj) a z , where a y denotes 
the annihilation operator of the mode of the flexural mo- 
tion in y-direction and the coupling strength is 

g' = <f> f I p L5 y /$o, (19) 
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where 8 y is the zero point fluctuation of the NAMR in 
the y-direction. However, Since 
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B 



(20) 



Bbias = &f/S, with S the area enclosed by the loop 
of 3- J J flux qubit, this additional coupling can be sub- 
stantially suppressed by a properly designed asymmet- 
ric structure of the NAMR to set 8 y /S z -C 1. This can 
be made when the dimension of the NAMR in the y- 
direction is larger than that of the z-direction. Then the 
dominant coupling is the one induced by the magnetic 
field B . The bias field for the flux qubit at ®f_ = 0.5<S> 
is calculated as 19 /xT with the l oop area of Ref [60] or 250 
//T with a smaller loop of Ref [40]. This corresponds a 
coupling strength of 12 kHz or 160 kHz (with S y = 0.15 Z ), 
which is always less than 10" 1 of the coupling induced 
by the external applied magnetic field Bq even at the 
"weak coupling" situation discussed above. Therefore, 
this always-on coupling is negligible for Bq stronger than 
1 mT. 

Another possible difficulty in experiments might come 
from the vibration of the sample (the substrate with the 
flux qubit and the NAMR on it). If the controlling mag- 
netic field -Bo in the y-direction is generated by a coil lo- 
cated on another cold finger, then in general, there may 
exist uncontrollable relative motion of this sample cavity 
against the outside coil system. This relative oscillation 
between the coil and the sample induces the fluctuation 
of the controlling field. And more seriously, the torsional 
oscillation of the sample will cause the deviation of the 
controlling magnetic field Bq away from y-direction, i.e., 
the angle 9 between Bq and the plane of the qubit-NAMR 
loop can not be zero strictly. Thus the none-zero com- 
ponent of Bq in z-direction BqSuyO ~ Bq6 penetrating 
the loop causes qubit energy fluctuations and decoher- 
ence. Roughly speaking, this decoherence source is pro- 
portional to Bq. Therefore, one need to optimize Bq such 
that it is large enough to dominate the coupling but not 
too large to induce strong decoherence from the torsional 
vibration of the sample plane with respect to Bq. 

On the other hand, we could also reduce the fluctua- 
tion of the bias magnetic field of flux qubit by improving 
the experimental setup. For example, we could use an 
8-shaped gradiometcr qubit. By the modified structure 
of the flux (as shown in figO, we can cancel the uni- 
form magnetic field fluctuations over the qubit in the 
z-direction. The magnetic field threading the loop of the 
flux qubit is measured by the readout SQUID. Then with 
suitable feedback to compensate the fluctuation men- 
tioned above, the bias magnetic field in the qubit-NAMR 
loop can be stabilized. In this case, strong magnetic field 
is attainable by the off-chip coil. However, scaling up to 
many qubits is not so straightforward and needs further 
consideration. 

Another possible solution to overcome the above ob- 
stacles is to prepare Bq by a superconducting coil and fix 
it on the sample chip at the dilution temperature. This 



B r 



B, 
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FIG. 7: (Color on line) The unexpected vibration induced 
fluctuation in the bias magnetic field could be canceled by 
feeding back the readout of the SQUID . 



method has advantageous that the setup of the proposed 
controllable coupling mechanism need not be modified 
and sufficient strong controlling magnetic field Bq can be 
achieved. 



VIII. DISCUSSIONS AND REMARKS 

In summary, we propose a novel solid-state cavity QED 
architecture that can reach the "strong coupling regime" 
based on a superconducting flux qubit and an NAMR. 
In this composite system, the quantized flexural mode of 
the NAMR is coupled with the persistent current gener- 
ated in the superconducting loop. The coupling strength 
can be independently modulated by an external magnetic 
field. We study the entangled energy spectrum of this 
composite system and find that a QND measurement of 
flux qubit can be made in the dispersive limit. This com- 
posite system can be scaled up and the coupling mecha- 
nism can be extended to the case that the superconduct- 
ing junction is a multi-level system. We also carefully 
examine the practical issues for experimental realization. 
The controllable strong coupling and the scalability en- 
able coherent control over this system for the quantum 
information processing as well as quantum state engi- 
neering. Besides, this cavity QED architecture offers a 
new scenario to demonstrate the intriguing phenomena 
of quantum optics in solid-state quantum device. It also 
provides a possibility to test the quantization effect of 
mechanical motion. 
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